Nonuniform Rashba-Dresselhaus spin precession along arbitrary paths 
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Electron spin precession in nonuniform Rashba-Dresselhaus two-dimensional electron systems along arbi- 
trary continuous paths is investigated. We derive an analytical formula to describe the spin vectors (expectation 
values of the injected spin) in such conditions using a contour-integral method. The obtained formalism is ca- 
pable of dealing with the nonuniformity of the Rashba spin-orbit field due to the inherent random distribution 
of the ionized dopants, and can be applied to curved one-dimensional quantum wires. Interesting examples are 
given, and the modification to the spin precession pattern in a Rashba-Dresselhaus channel when taking the 
random Rashba field into account is shown. 
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I. INTRODUCTION 

Space-inversion-asymmetry-induced spin precession in 
two-dimensional electron gas (2DEG) channels opens the pos- 
sibility of spin-field-effect transistor (SFET) (Ref. 1) and is it- 
self an intriguing physical phenomenon in the emerging field 
of semiconductor spintronics.^"' The heart of such spin preces- 
sion lies in the well-known Kramers theorem,'* which predicts 
an intrinsic zero-field spin-splitting in inversion-asymmetric 
crystal structures in the absence of applied magnetic fields. 
The corresponding spin-splitting can be described by an ef- 
fective magnetic field, forcing the transport spins to precess 
about such fields. 

Mechanisms of inversion asymmetry in 2DEG systems 
include structure inversion asymmetry, more commonly re- 
ferred to as the Rashba spin-orbit (RSO) interaction,^ of the 
confining potential, and the bulk inversion asymmetry, re- 
ferred to as the Dresselhaus spin-orbit (DSO) term for zinc- 
blende semiconductors,^ of the underlying crystal structure. 
Specifically, the Rashba term originates from the transforma- 
tion of the electric field felt by the electron to an effective 
momentum-dependent magnetic field. Hence factors influenc- 
ing the RSO strength involve local charge distribution (such as 
charged impurities^ or ionized dopant atoms** ''), and applied 
external electric fields (usually the gate voltage). 

Perhaps the gate-voltage tunability of the RSO strength, 
quite successfully demonstrated in InGaAs/InAlAs-based 
(Ref. [13) and GaAs/AlAs-based (Refs. [O and [H) quan- 
tum wells (QWs), determines the precession rate of the in- 
jected electron spin, and is therefore the most important key 
to the realization of the SFET. However, more than a decade 
of struggle has still yielded disappointing results. Its failure 
may be ascribed to the so far less successful factors:^ (i) ef- 
ficient spin injection rate and (ii) uniformity of the RSO in- 
teraction. Apart from the most challenging issue — the spin 
injection problem, nonuniformity of the RSO coupling aris- 
ing from the inherent random distribution of ionized dopants 
in the vicinity of the 2DEG plane has recently attracted cer- 
tain attention to reexamine the Rashba 2DEG systems.^^iS This 
implies that the RSO coupling strength a should be treated 
as position dependent, and the corresponding modification to, 
e.g., spin precession rate, in the Rashba channel must be taken 
into account. 



In this paper we construct a contour-integral method to deal 
with spin precession in nonuniform Rashba-Dresselhaus (RD) 
2DEGs along arbitrary paths. We derive an analytical formula 
to describe the spin vectors, i.e., the expectation values of spin 
operators done with respect to the injected electron spin state, 
inside the 2DEG channel. The obtained formula is found to 
be a direct generalization of the previous one,'^ which is ap- 
plicable to injected spins with straight space evolution in uni- 
form RD 2DEGs. The constructed formalism can straight- 
forwardly deal with the nonuniformity of the RSO coupling, 
DSO coupling, and even the effective mass m*, and is appli- 
cable to curved one-dimensional (ID) channels, e.g., quantum 
ring structures. Throughout this paper, single-particle quan- 
tum mechanics is used and the effective mass approximation 
is adopted. 



II. FORMULA 

We begin by considering a [001] -grown zinc-blende-based 
2DEG, where the Rashba and Dresselhaus Hamiltonians 
are described by Hn = a/h{pycr^ ~ pxcr^) and Hd = 
j3/h {pxcr^ — PyO'^), respectively. Whereas the RSO strength 
a stems from the electric field generated by the ionized 
dopants in the vicinity of the 2DEG plane, a is in principle 
position dependent and material specific. On the other hand, 
the DSO strength, in bulk systems, arises from the inversion- 
asymmetry of the crystal structure and thus may not vary 
with different positions if the sample is well grown. How- 
ever, when restricted to 2DEGs, the DSO parameter f3 de- 
pends on the width of the confining quantum well, and may 
therefore also be position dependent. When the underlying 
crystal structure of the 2DEG channel is not perfectly grown, 
spatially varying effective mass (r) may also occur. 

In the uniform case, where all of the three parameters a, j3, 
and m* are constant in space, an electron spin injected at 
can be described by (see Appendix A) 



(1) 



when detected at after a straight space evolution ri Vd- 
Here, |s)- is a given spinor describing the spin injected at 
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FIG. 1: (Color online) Schematic of a continuous curved path sec- 
tioned into pieces. 



■Ti, kis the average of the two projected spin-dependent wave 
vector given by Eq. ( IA4I ). and 

|V'<T;0r-d-n) = ^ j (2) 

is the RD eigenspinor with 4>ra~ri being the argument 
of the displacement vector — and f (a, /3, 0) = 
arg [(a cos (/) + /3sin (/)) + i (a sin0 + /3 cos (/))]. The phase 
difference is defined by 

A0P = ^"^''C (-^r) r , (3) 

where 

C = + 2a/3sin(20f.) (4) 

is the composite spin-orbit coupling strength. Note that in this 
paper the spatially evolved state ket without superscript, say 
\s)p^^p^, means space evolution from position ri to position 
r2 through the straight path connecting the two points, while 
|s)? means evolution through path C. 

For the nonuniform case, we consider an arbitrary continu- 
ous path C connecting the injection and detection points, and 
section the path into N pieces (see Fig. [T]i. In the limit of 
N ^ oo, each section approaches to a straight path with con- 
stant a, (3, and m*, and Eqs. ([TJ and (|3]l are then applicable 



for each section. Again, we inject a spin at and detect the 
spin at frf through, however, the curved path C along which 
the parameters a, f3, and m* may be functions of position. Af- 
ter successive application of Eq. ([T]|, we obtain (see Appendix 
B) 

where the total phase difference is given by a contour integral 

AQ = ^ [ m*Cdf. (6) 
" Jc 

Compared to the uniform case Eqs. ([B and (|3]l, the problem of 
nonuniform systems is simply extended to integrate the total 
phase difference encountered by the electron along the path 
C. 

We remind here that the above formalism (and later the de- 
rived spin vector formula) is applicable for nonuniform sys- 
tems but only with continuous spatial evolution. Propagation 
along discontinuous paths must be handled by sectioning the 
path into pieces of continuous ones [since Eq. ( IB6I ) does not 
hold on discontinuous points] . Another discontinuity that may 
crash the above formalism is the reverse of the RD field di- 
rection [at which Eq. ( IB6I ) also fails]. In pure Rashba (such 
as Si-Si2:Gei_2: asymmetric QWs) cases, the change of sign 
a ^ —a, due to either the random dopants or the reverse gat- 
ing induces such discontinuity. In the composite case (both 
a and /? nonvanishing), the cancellation of the momentum 
dependence of the RD effective magnetic field at the condi- 
tion \a\ ~ 1/3 1 provides possibility to run the SFET in the 
nonballistic regime'^ and also contains fundamental physical 
phenomena. However, the corresponding field direction may 
become unstable along ±[110] or ±[110] directions, even 
though the field magnitude is still continuous. When dealing 
with such systems, one must carefully check the continuity of 
the RD field direction and do the calculation partitively. 

Now we apply Eq. (|5]l to perform the spin expectation 
values (S) = h/2 (a), where a is the Pauli matrix vector 
[a^, ay^az)- After some mathematical manipulation, we ob- 
tain 



cos 6s cos ipd sin A0 + sin Og 
- cos 9a sin ipd sin + sin 9s 



^ N 2 ^0 , ^ ^ ■ 2 Ae" 

cos [ifd - ipi + (ps) COS — COS [Lfd + '■Pi - <Ps) sm ^- 

. ^ 2 Ae . , ^ , . 2 Ae 

sm [ipd - ipt + (ps) cos^ — sm (ipd + (fi - (ps) sm ^- 



(7) 



cos 9s cos A8 + sin 9s cos [ipi — (ps) sin AO 
I 



with ipi(^d) = arg[(a^(d) cos0i(d) + /3j(d) sin 4(d)) + 
j(a,(d)Sin0,(d) + /3,(d) cos 4(d))], 9s and cps being the po- 
lar and azimuthal angles of the injected spin, respectively. 



Clearly, Eq. (O recovers the previous spin vector formula— 
when a, /3, m*, and (p are constants, i.e., uniform spin preces- 
sion along straight paths. 
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III. RESULTS AND DISCUSSION 

In this section we apply our generalized version of spin- 
vector formula Eq. dTji with the total phase AQ given by Eq. 
^ first to some interesting ID geometries, including quan- 
tum rings and quantum wires, and then to a more realistic 
2DEG case, where position dependent RSO coupling, gener- 
ated by regular and random distributions of ionized dopants, 
is present. In all cases, we assume constant DSO coupling and 
the electron effective mass. 
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A. ID quantum ring 

Beginning with an ideal Rashba half-ring with ring radius 
0.1 /im, we inject an electron at one end, and detect (S) down 
the way along the ring to the other end. The half-ring is as- 
sumed to be ideally ID and made of InGaAs-based materials. 
We set the RSO parameter a — 0.03 cVnm and the electron 
effective mass m* = 0.03mei^ Two configurations of the 
injected spin are considered. In ring (i) of Fig. Ha), we in- 
ject the spin with polarization parallel to the radial direction, 
which is also the eigenstate of this Rashba ring case. As ex- 
pected, the spin vectors maintain in the radial direction since 
the spin is injected in its eigenstate and will stay in this state 
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FIG. 2: (Color online) (a) Spin vectors along four ideal ID half- 
rings. Rings (i) and (ii) contains only the RSO while in rings (iii) 
and (iv) the DSO is involved. The spin, polarized parallel either to 
the radial or the tangential directions, is injected at the lower end of 
each ring, (b) Corresponding (Sz) in units of ft/2 as a function of 
the ring argument 9. 



FIG. 3: (Color online) Two ideal ID wires with position-dependent 
RSO strength. The upper wire shows the effect of a linear a and the 
lower wire shows a sine-varying a. In the main panel, a and (Sz) as 
functions of the longitudinal position are plotted. 



hereafter. This can also be viewed as a generalized version of 
our previous proposal of precessionless spin transport wire. — 
In Rashba systems, one can precessionlessly transport spins 
even through an arbitrary wire shape (but with continuous cur- 
vature), once the spin is injected with polarization perpendic- 
ular to the wire direction. 

Next we inject the spin parallel to the tangential direction. 
As shown in ring (ii) of Fig. 12 a), the injected spin precesses 
upright down the way to the end. This is equivalent to the 
Datta-Das SFET with a curved ID Rashba channel. Whereas 
the RSO coupling induces an effective magnetic field i?eff, 
which is always perpendicular to the electron transport (Sett | 
the radial direction here), rings (i) and (ii) in Fig. |2a) are both 
reasonable and expectable. However, when the DSO coupling 
is involved, the z-rotational symmetry is broken, and the spin 
vectors cannot be expected intuitively, but can still be well de- 
scribed by our formula Eq. Q [see rings (iii) and (iv) in Fig. 
|2ja)]. To clarify the influence due to the DSO, we plot (Sz) 
as a function of the ring argument in Fig. |2jb), where line 
(i) and line (ii) exhibit precessionless and complete precess- 
ing behaviors, respectively, and change to line (iii) and line 
(iv) after getting the DSO coupling involved. 

B. ID quantum wire 

We now consider two straight quantum wires with position- 
dependent RSO strength, one being a linear a wire and the 
other a sine-varying a wire. These wires exhibit interesting 
spin precession behaviors, as shown in Fig. [3] where the up- 
per and lower wires, sketched from the side view, correspond 
to the Unear and sine-varying Rashba wires, respectively. As 
expected, the injected spin rotates with a precession rate in- 
creasing with the position, in the linear Rashba wire. In the 
sine case, the alternate positive and negative RSO drives the 
spin to precess backward (when a > 0) and forward (when 
a < 0). These quantum wires may be fabricated by spe- 
cial gating designs, but applicability to the spintronic device 
is not obvious. However, the physics indicated here is sim- 
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pie and clear. Whereas the magnitude of the RSO strength a 
determines the spin precession rate, the sign of a controls the 
precession direction. From the viewpoint of gating, a strong 
enough backgate may reverse the electric field direction in the 
2DEG, and hence the resulting Rashba field direction, forcing 
the spin to precess in an opposite direction. From the view- 
point of the random Rashba field, such an alternate a is similar 
to that in a symmetric quantum well, e.g., a Si-Sia;Gei_a; well, 
and will induce a finite spin relaxation ratei^ 



C. 2DEG channel with regular and random dopants: Random 
Rashba effect 

Finally, we arrive at the crucial but realistic problem: how 
will the spin precession pattern be modified due to the random 
Rashba field? A first guess on the basis of Eq. Q is that how 
wide the angle the spin precesses through lies in how much 
spin-orbit interaction the electron encounters along the path it 
goes. A symmetric well, which has a vanishing spatially av- 
eraged RSO (a) — 0, does not have spin precession pattern 
within a submicrometer scale, even though there are nonvan- 
ishing fluctuating Rashba fields. This is similar to the sine- 
varying Rashba wire shown in Fig. [3]but with a much rapider 
a oscillation period atomic sizes). We now consider an 
asymmetric well, and adopt the same principle introduced by 
Sherman et al.'^ to generate the random Rashba field. 

Specifically, we consider a 400 nm x 267 nm InAs-based 
2DEG with only a single dopant layer, located at a spacing 
zo — 20 nm from the conducting plane. Within the dopant 
layer, charged dopant atoms are randomly distributed, gen- 
erating the Coulomb field with z component described by 



100 200 300 I 



-3/2 



Here e is the 



electron charge, e is the material permittivity, p = {x, y) 
is the two-dimensional position vector, and fj is the posi- 
tion of the jth dopant. We take the dopant concentration 
n = 2.5 X 10^^ cm^-^ (same with Ref. the permittiv- 
ity e = 15.1eo (eo is the free space permittivity), and the 
lattice constant a — 0.6058 nm.'^ Using the linear model 
a{p) = asocEz [p) with aso = 110 for InAs,!iii2 we 
obtain the Rashba field in the cases of regular and random 
dopant distributions in Figs. Ufa) and|4fb), respectively. 

In the 2DEG channel, we also take into account the DSO 
coupling by assuming the well thickness w = 50 A, leading 
to (3 ^l{kl) « 1.062 X 10-3 eVnm. Here 7 = 26.9 eVA^ 

for InAs QWs (Ref. 11) and (fc^) « {n/w)'^ in the case of 
rigid walls. Assuming point injection of a; -polarized spins at 
the left-center of the channel with transport direction [110], 
the spin vector formula Eq. d?) gives the spin precession pat- 
terns, as shown in Fig. |4|c) for the regular channel and Fig. 
Hid) for the random channel. Clearly, the randomness of the 
dopant distribution distorts the precession pattern, and hence 
is expected to lower the purity of the spin signal collected 
at the drain end. This may bring another difficulty in real- 
izing the Datta-Das transistor. ' Moreover, such unpredictable 
random distribution will eventually cause an uncertainty of 
the collected spin signal, for each individual 2DEG channel. 
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FIG. 4: (Color online) Position-dependent Rashba field a with (a) 
regular and (b) random dopant distribution, resulting in the corre- 
sponding spin precession patterns in (c) and (d), respectively, where 
point injection of a; -polarized spins at the left ends, are assumed. The 
color bars determine the magnitude of a in units of eV nm in (a) and 
(b), and the z compont of the spin vector (Sz) in units of h/2 in (c) 
and (d). (e) and (f) show (Sx) , {Sy) , (Sz) , and a as functions of 
longitudinal position x at y = for the regular and the random cases, 
respectively. 



even if the dopant concentration is perfectly controlled. For 
example, the standard deviation calculation with the 2DEG 
channel conditions considered here shows that AS^:/ {fi/2) ~ 
O (lO-i) for InAs-based and AS^/ {h/2) ~ O (lO^^) for 
GaAs-based materials. (Here ASi means the standard devia- 
tion of the averaged i component of the spin signal collected 
at the drain end in the channel, among a number of different 
samplings.) 

Note that we have chosen [110] as the channel direction, 
where the Rashba and Dresselhaus fields are parallel,'^ and 
the random effect seems moderate. When analyzing other 
channel directions such as [110], the obtained spin precession 
pattern will be totally different from the regular case. Last, if 
we go back to the ID case and focus on the center path from 
the source end to the drain end, the distortion of the spin pre- 
cession curves due to the random Rashba field is surprisingly 
weak [see Figs. |4je) andUf)]. This again supports the sug- 
gestion of using ID or quasi- ID channels for the Datta-Das 
transistor. 



IV. SUMMARY 

In conclusion, we have derived an analytical formula to de- 
scribe the spin precession in nonuniform RD 2DEGs, using a 
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contour-integral method. The obtained results are applicable 
for 2DEG systems with position-dependent a, (3, and to*, and 
curved ID wires. Numerical examples show interesting spin 
precession behaviors, and also support the idea of the preces- 
sionless spin transport wire'^ in the Rashba ring case. We 
have also demonstrated how the spin precession pattern in a 
RD 2DEG may be distorted when taking the random Rashba 
field into account. In the viewpoint of carrying out the SFET, 
the random Rashba effect may distort the regular spin pre- 
cession pattern in the 2DEG channel, especially along weak 
spin-orbit strength directions. Such undesired influence can 
be suppressed quite well when narrowing down the channel to 
ID [compare Figs. lUf) with|4fe)], supporting the pioneering 
suggestions of using quasi- ID channels.''^" 



Acknowledgments 

One of the authors (M.H.L.) is grateful to Son-Hsien Chen 
and Ivo Klik for valuable discussions and suggestions. This 
work is supported by the Republic of China National Science 
Council Grant No. 94-2112-M-002-004. 



APPENDIX A: UNIFORM SPACE EVOLUTION ALONG 
STRAIGHT PATH 

To derive Eq. ([T), we use the translation operator to bring 



the ket Is)-, from ri to rd- 



Jp-{rd-ri)/h 



Expanding |s) - in terms of the RD eigenspinors \iIj±; 
and denoting Ar = jr^ — ri |, we have 

= e^P^^/'^ (c+ |^+; + c_ IV^-; , 

(Al) 

where c± — {tp±; 0f^_fjs)- are the expansion coefficients 
(0f is the argument angle of the vector r). Note that the elec- 
tron carrying momentum p is supposed to move in the di- 
rection rd — Ti, so the product p ■ {fd — ri) is simplified to 
p\rd — ri\ = pAr. Defining 



- k+ = Ak = 



2m*( 



(A2) 



where C, is the composite spin-orbit strength given in Eq. (HI, 
we rewrite Eq. (lAlb as 



Ak- Ar 



APPENDIX B: NONUNIFORM SPACE EVOLUTION 
ALONG CURVED PATH 

To derive Eqs. dD and (|6]l, we refer to Fig. [T| where we 
denote ri(^d) = ^^1(^+1) successively apply Eq. ([T]i from 
the last to the first section. Taking the jxh section for example, 
we have 



(Bl) 



cr=±l 



with Ar = 



Ae, ^ 



2m*Cj^r/h^, and Q = \J o?j + /3| + 2ajPj sin (2(/)j). Note 
that the RSO strength, DSO strength, and effective mass are 
ideally local: aj = a (fj), f3j = (3 (r*,), and m* = m* (rj). 
Note also that the common phase factor outside the summa- 
tion exp (ifcj A?-) with kj = {^k^ + ^7) depend not 
only on the position j^j but also on the direction (pj due to the 
anisotropic Fermi contour in the general RD case. However, 
we will show that, as in the uniform straight case, this phase 
does not contribute to the spin vector (S). In the following 
we will abbreviate the RD eigenspinor lipa-', 4'j) to lipt) for 
simplicity. 

For the fuU translation, we can start with the last section 



a=±l 



(B2) 



and successively substitute 



|„\ _ ifeiv-iAr 



cr=±l 

/,A'-1|„\ L;,JV-1\ 

' nv-2^rjv- 



I A _ ifejv-2Ar 

P/fjv-2^rjv_l 



E 



i(TAejv_2/2 



cr=±l 



(B3a) 



(B3b) 



cr=±l 



(B3c) 



_ gi/cAr y^ g-i(TA0/2 
cr=± 



with 



Cct IV'ct; ^Fd-fi) J (A3) into the expression Eq. (IB2I) . Then totally we have 

E 



cr(«).cr(«-l). 



k = 



(A4) 



and A^^ = AfcAr. Note that the common phase e^'"^^ will 
vanish when calculating expectation values. Whereas k de- 
pends on the energy the electron spin carries while Afc, or the 
phase difference A6, depends only on the spin-orbit strength, 
the physics of |s)-^-^, such as the spatial spin precession, 
does not depend on the injection energy. 



^ g-i<T(")Aejv/2g-i<T("-i'Aejv_i/2 , , , g-j<T"Ae2/2 

X |V<x') • • • |V'^(iV-2)) |V^^(iV_l,)) |V'<,(iV)) , (B4) 



where we have further abbreviated 



to 



|?/'cr(3))- One must bear in mind that the short- 
hand lipcrii)) actually means \tp^(j); (f>j). In Eq. 
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(IB4l i. the spinor overlaps can be rearranged as 
For the jth bracket, the overlap reads 

(V.u) IV-.o-d) = ^ (e^(^^-^-i) + (T«a(^-i)) . (B5) 

In the limit of iV ^ cx3 and the continuous condition, 
this overlap becomes 5^{j) ^(j -1) since ~ (Xj ~ 

ctj-i, and /3j ~ [3j^i yield (^^ w '■Pj-i^ so that we have 
« (H-(T(^)a(J-i))/2 = ,5,0)_,„-i). There- 
fore, the spinor overlaps are reduced to 

('0<T'k)o (■0<T"|V'(T') ■ • • ('0^(iV-l)|V'<x(«-2)) (7/'^(iV)|7/'^(iV-l)) 

= '^(t("),(t(«-i)'^(t(«-i),ct("-2) • ■ ■ (^(t",!!' (V'ct' |'S)o ' (1^6) 

leading to a greatly simplified expression 



exp j i A:„ Ar j 

\ n=l / 



N 



(B7) 



where we have returned the notations for the RD eigenspinors 
= IV'cr') and \ti}a;4>d) = iV'o-*"))' and also the input 
spin |s)-, = |s) Rewriting the prefactor and the total phase 



into integration forms 

expi|^^fc„Arj ^ expi / fc (</>) dr , (B8) 

\n=l / p 

^ ^ 2m*C- 2 /■ 

E^^^=-E-?^^''^:a^/"^*Crfr^Ae, (B9) 

we finally obtain 

|s)? - = e*-''c ^("^)'^'' 



(BIO) 

Since the prefactor does not survive when doing expectation 
values (such as (it)) with the above ket, it makes no difference 
to express |s)?^-^ as Eq. (|5]). 

Although we have shown that the common phases in Eqs. 
(IA3b and ( IB 1 Oi l are irrelevant when calculating, e.g., the spin- 
vectors {S), we remind here that these phases will become 
relevant once superposition of different states is involved. In 
that case the prefactor induces oscillations due to interference 
between the superposed states, and is therefore not negligible, 
unless the inverse Fermi momentum is much larger than the 
system size, i.e., fcAr <C 1. 
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